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Abstract We consider the modification of the formulas for black hole radiation, due 
to the self-gravitation of the radiation. This is done by truncating the coupled particle- 
hole system to a small set of modes, that are plausibly the most significant ones, and 
quantizing the reduced system. In this way we find that the particles no longer move 
along geodesies, nor is the action along the rays zero for a massless particle. The radiation 
is no longer thermal, but is corrected in a definite way that we calculate. Our methods 
can be extended in a straightforward manner to discuss correlations in the radiation, or 
between incoming particles and the radiation. 



1. Introduction 

Black hole radiance |]J was originally derived in an approximation where the back- 
ground geometry was given, by calculating the response of quantum fields to this (collapse) 
geometry. In this approximation the radiation is thermal, and much has been made both 
of the supposed depth of this result and of the paradoxes that ensue if it is taken literally. 
For if the radiation is accurately thermal there is no connection between what went into 
the hole and what comes out, a possibility which is difficult to reconcile with unitary 
evolution in quantum theory - or, more simply, with the idea that there are equations 
uniquely connecting the past with the future. To address such questions convincingly, one 
must go beyond the approximation of treating the geometry as given, and treat it too as 
a quantum variable. This is not easy, and as far as we know no concrete correction to the 
original result has previously been derived in spite of much effort over more than twenty 
years. Here we shall calculate what is plausibly the leading correction to the emission 
rate of single particles in the limit of large Schwarzschild holes, by a method that can be 
generalized in several directions, as we shall outline. 

There is a semi-trivial fact about the classic results for black hole radiation, that 
clearly prevents the radiation from being accurately thermal. This is the effect that 
the temperature of the hole depends upon its mass, so that in calculating the "thermal" 
emission rate one must know what mass of the black hole to use - but the mass is different, 
before and after the radiation! (Note that a rigorous identification of the temperature of 
a hot body from its radiation, can only be made for sufficiently high frequencies, such 
that the gray-body factors approach unity. But it is just in this limit that the ambiguity 
mentioned above is most serious.) As has been emphasized elsewhere 0, this problem 
is particularly quantitatively acute for near-extremal holes — it is a general problem for 
bodies with finite heat capacity, and in the near-extremal limit the heat capacity of the 
black hole vanishes. 

To resolve the above-mentioned ambiguity, one clearly must allow the geometry to 
fluctuate, namely to support black holes of different mass. Another point of view is that 
one must take into account the self- gravitational interaction of the radiation. 



2. Model and Strategy of Calculation 

To obtain a complete description of a self-gravitating particle it would be necessary to 
compute the action for an arbitrary motion of the particle and gravitational field. While 
writing down a formal expression for such an object is straightforward, it is of little use 
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in solving a concrete problem due to the large number of degrees of freedom present. 
To arrive at a more workable description of the particle-hole system, we will keep only 
those degrees of freedom which are most relevant to the problem of particle emission from 
regions of low curvature. The first important restriction is made by considering only 
spherically symmetric field configurations, and treating the particle as a spherical shell. 
This is an interesting case since black hole radiation into a scalar field occurs primarily 
in the s-wave, and virtual transitions to higher partial wave configurations are formally 
suppressed by powers of t^. 

Before launching into the detailed calculation, which becomes rather intricate, it seems 
appropriate briefly to describe its underlying logic. After the truncation to s-wave, the 
remaining dynamics describes a shell of matter interacting with a black hole of fixed mass 
and with itself. (The mass as seen from infinity is the total mass, including that from 
the shell variable, and is allowed to vary. One could equally well have chosen the total 
mass constant, and allowed the hole mass to vary.) There is effectively one degree of 
freedom, corresponding to the position of the shell, but to isolate it one must choose 
appropriate variables and solve constraints, since the original action superficially appears 
to contain much more than this. Having done that, one obtains an effective action for the 
true degree of freedom. This effective action is nonlocal, and its full quantization would 
require one to resolve factor-ordering ambiguities, which appears very difficult. Hence 
we quantize it semi-classically, essentially by using the WKB approximation. After doing 
this one arrives at a non-linear first order partial differential equation for the phase of the 
wave function. This differential equation may be solved by the method of characteristics. 
According to this method, one solves for the characteristics, specifies the function to be 
determined along a generic initial surface (intersecting the characteristics transversally) , 
and evolves the function away from the initial surface, by integrating the action along the 
characteristics. (For a nice brief account of this, see 0.) 

When the background geometry is regarded as fixed the characteristics for particle 
motion are simply the geodesies in that geometry, and they are essentially independent 
of the particle's mass or energy — principle of equivalence — except that null geodesies 
are used for massless particles, and timelike geodesies for massive particles. Here we 
find that the characteristics depend on the mass and energy in a highly non-trivial way. 
Also the action along the characteristics, which would be zero for a massless particle 
and proportional to the length for a massive particle, is now a much more complicated 
expression. Nevertheless we can solve the equations, to obtain the proper modes for our 
problem. 

Having obtained the modes, the final step is to identify the state of the quantum field 

* Since we do not address the ultraviolet problems of quantum gravity these corrections are actually 
infinite, but one might anticipate that in gravity theory with satisfactory ultraviolet behavior the virtual 
transitions will supply additive corrections of order w 2 /A 2 , where A is the effective cutoff, but will not 
alter the exponential factors we compute. 
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- that is, the occupation of the modes — appropriate to the physical conditions we wish 
to describe. We do this by demanding that a freely falling observer passing through the 
horizon see no singular behavior, and that positive frequency modes are unoccupied in 
the distant past. This, it has been argued, is plausibly the appropriate prescription for 
the state of the quantum field excited by collapse of matter into a black hole, at least in 
so far as it leads to late-time radiation. Using it, we obtain a mixture of positive- and 
negative- frequency modes at late times, which can be interpreted as a state of radiation 
from the hole. For massless scalar particles, we carry the explicit calculation far enough 
to identify the leading correction to the exponential dependence of the radiation intensity 
on frequency. 



3. Effective Action 

We now derive the Hamiltonian effective action for a self-gravitating particle in the 
s-wave. The Hamiltonian formulation of spherically symmetric gravity is known as the 
BCMN model [|[]; our treatment of this model follows that of ||. First, we would like 
to explain why the Hamiltonian form of the action is particularly well suited to our 
problem. As explained above, our physical problem really contains just one degree of 
freedom, but the original action appears to contain several. The reason of course is that 
Einstein gravity is a theory with constraints and one should only include a subset of the 
spherically symmetric configurations in the physical description, namely those satisfying 
the constraints. In general, in eliminating constraints Hamiltonian methods are more 
flexible than Lagrangian methods. This appears to be very much the case for our problem, 
as we now discuss. 

In terms of the variables appearing in the Lagrangian description, the constraints have 
the form 



r, r; 



0. 



where f is the shell radius, and ' represents 4. When applied to the spherically symmetric, 
source free, solutions, one obtains the content of Birkhoff 's theorem - the unique solution 
is the Schwarzschild geometry with some mass, M. Since this must hold for the regions 
interior and exterior to the shell (with a different mass M for each), and since M must be 
time independent, we see that only those shell trajectories which are "energy conserving" 
are compatible with the constraints. This feature makes the transition to the quantum 
theory rather difficult, as one desires an expression for the action valid for an arbitrary shell 
trajectory. This defect is remedied in the Hamiltonian formulation, where the constraints 
are expressed in terms of momenta rather than time derivatives, 

Ch [r,p;gij,7Tij] = 0. 
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At each time, the unique solution is again some slice of the Schwarzschild geometry, but 
the constraints no longer prevent M from being time dependent. Thus, an arbitrary shell 
trajectory f(i), p(t), is perfectly consistent with the Hamiltonian form of the constraints, 
making quantization much more convenient. 

The starting point for the Hamiltonian formulation of gravity is to write the metric 
in ADM form gf: 

ds 2 = -N l {t, rfdt 2 + L(t, rf[dr + N r {t, r)dtf + R(t, r) 2 [d6 2 + sinfl 2 d0 2 ] (3.1) 

In considering the above form, we have restricted ourselves to spherically symmetric 
geometries at the outset. With this choice of variables, the action for the shell is written 

S s = -m J ^-g^dx^dx" = -mjdt \Jn i2 - L 2 (r + JV r ) 2 , (3.2) 

m representing the rest mass of the shell, and the carets instructing one to evaluate 
quantities at the shell (jj^ = ^(i, f)J. 

The action for the gravity-shell system is then 
S = — [d 4 x y/^gTZ -m fdt \l (N 1 ) 2 - L 2 (f + N r ) 2 + boundary terms (3.3) 

167T J J 

and can be written in canonical form as 
S = Jdtp'r + Jdtdr [vr R R + vr L L - N\H s t + Hf) - N r {H s r + H?)} - Jdt M ADM (3.4) 

with 

—p5{r — r) (3.5) 

U G r = R'tt r - Ln' L (3.6) 

where ' represents 4-. Madm is the ADM mass of the system. The inclusion of this 
last term deserves some comment. Because the Einstein-Hilbert action contains second 
derivatives, a general variation of the metric variables gives a nonzero result even when the 
equations of motion are satisfied and the variation of the metric is zero on the boundary 
of the space. If we restrict the class of metrics to those which are asymptotically flat, 
with A^* — > 1 and A^ r — > as r — ► oo, then the variation of the term / dtM^DU. precisely 
cancels the unwanted terms, and so gives a well defined variational principle 0. It is 
important to note that Madm is a function of the metric variables (whose explicit form 
will be displayed in due course) and is numerically equal to the total mass of the combined 
gravity-shell system. 



U\ = \l{p/L) 2 + m 2 8{r -f) ■ H s r = 

Ln L 2 n L n R ( RR'Y R' 2 L 
2R 2 R \ L ) 2L 2 ' 
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We now wish to eliminate the gravitational degrees of freedom in order to obtain an 
effective action which depends only on the shell variables. To accomplish this, we first 
identify the constraints which are obtained by varying with respect to N* and N r : 

Ht = H\ + H? = ; H r = H s r + = 0. (3.7) 

By solving these constraints, and inserting the solutions back into (£1) we can eliminate 
the dependence on n R and til- We first consider the linear combination of constraints 

= %-Ut + ^-H r = —M' + S-Ht + ^H s r (3.8) 
L RL L RL 

where 

M = ^ + *-^. (3.9) 
2R 2 2L 2 v ' 

Away from the shell the solution of this constraint is simply A4 = constant. By considering 

a static slice (ttl = ttr = 0), we see that the solution is a static slice of the Schwarzschild 

geometry with Ai the corresponding mass parameter. The presence of the shell causes 

M. to be discontinuous at f, so we write 

M = M r < r 

M = M+ r>r. (3.10) 



As there is no matter outside the shell we also have Madm = M + . Then, using (|3.8|) and 
Q3.9Q we can solve the constraints to find til and n R : 

n L = RyJ{R'/L)* - 1 + 2M/R ; 7t R = ^ir' L r<r 

L 



tt l = R^{R'/Lf - 1 + 2M + /R ; n R = —7r' L] r > f. (3.11) 

The relation between M + and M is found by solving the constraints at the position of the 
shell. This is done most easily by choosing coordinates such that L and R are continuous 
as one crosses the shell, and are free of singularities there. Then, integration of the 
constraints across the shell yields 

?T£ if" + e) - 7r L (f - e) = -p/L 



R'(f + e) - R'(f - e) = -^\]p 2 + m 2 L 2 (3.12) 

R 

Now, when the constraints are satisfied a variation of the action takes the form 

dS = pdr + J dr(7r R 5R + ix L 5L) - M+ dt (3.13) 
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where 7Tr } l are now understood to be given by (|3.11|) , and M + is determined by solving 
( 3.12|) . We wish to integrate the expression (|3.13| ) to find the action for an arbitrary shell 
trajectory. As discussed above, the geometry inside the shell is taken to be fixed (namely, 
M is held constant) while the geometry outside the shell will vary in order to satisfy the 
constraints. It is easiest to integrate the action by initially varying the geometry away 
from the shell. We first consider starting from an arbitrary geometry and varying L until 
ti\r = ttl = 0, while holding r,p, R, L fixed: 

JdS = J™ dr tf =0 5L n L 



I r r 2 dr R\J (R'/L) 2 - 1 + 2M/R + J™ € dr tf =0 6L RyJ (R'/L) 2 - 1 + 2M + /R 

'dr 



RLJ(R'/L) 2 - 1 + 2M/R + RR' log 



R' /L-y/(R' /L) 2 -1+2M/R 



^\1-2M/R\ 



~°° dr 



RLJ{R'/Lf 



2M+/R +RR' log 



R'/L- sj (R'/L) 2 -1+2M+/R 

y/\l-2M+/R\ 



where the lower limit of integration, 



(3.14) 

properly extends to the collapsing matter 
forming the black hole; its precise value will not be important. We have discarded the 
constant arising from the lower limit of the L integration. In the next stage we can vary 
L and i?,while keeping tir^l = 0, to some set geometry. Since the momenta vanish, there 
is no contribution to the action from this variation. 

It remains to consider nonzero variations at the shell. If an arbitrary variation of L 
and R is inserted into the final expression of ( |3.14| ) one finds 



dS = dr [tt r 5R + n L 8L] — 



dS 
OR' 



f + e - 



dS 
OR' 



r - e 



dS 

dR+—dM + . (3.15) 



Since R' is discontinuous at the shell, 

dS . . dS . 

— — (r + e) — (r — e) 

dR' dR' 

is nonvanishing and needs to be subtracted in order that the relations 



6S_ 

Jr 



5S_ 
6L 



will hold. From ( |3.14 ), the term to be subtracted is 
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~dS 
OR' 



r + e 



dS 
OR' 



(f-e) 



dR 



-dRRlog 



+ dRRlog 



R'if - e)/L - ^(R'(f-e)/L) 2 -l + 2M/R 
1 - 2M/R\ 



R'{f + e)/L - y/(R'(f + e)/L) 2 - 1 + 2M/R 



(3.16) 

\\-2M+/R\ 

Similarly, arbitrary variations of L and R induce a variation of M + causing the appearance 
of the final term in ( |3.15| ) . Thus we need to subtract 



dS 



■dM, 



* ^{R/LY-l + 2M + /R 

ar L- ; dM . 

1 - 2M + /R 4 



(3.17) 



Finally, we consider variations in p, f, and t. t variations simply give dS = —M + dt. We 
do not need to separately consider variations of p and f, since when the constraints are 
satisfied their variations are already accounted for in our expression for S, as will be 
shown. 

Collecting all of these terms, our final expression for the action reads 



S 



dr 



RLy (R'/L) 2 - 1 + 2M/R + log 



R'/L - J {R'/L) 2 - 1 + 2M/R 



1 - 2M/R\ 



+ / dr 



RLJ{R'/L) 2 - 1 + 2M+/R + RR' log 



R'/L - J (R'/L) 2 - 1 + 2M+/R 



1 - 2M + /R\ 



- dtRR 



log 



R'(f - e)/L - J(R'(f - e)/L) 2 - 1 + 2M/R 



1 - 2M/R\ 



+ log 



R'(f + e)/L - J(R'(f + e)/L) 2 - 1 + 2M/R 



1 - 2M + /R\ 



f roo LJ(R'/L) 2 - 1 + 2M+/R . r 

+ I"' L dr l-2M +/R M + " /* M - (3 ' 18 > 

To show that this is the correct expression we can differentiate it; then it can be seen 
explicitly that when the constraints are satisfied ( |3.13j ) holds. 
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We now wish to write the action in a more conventional form as the time integral of 



a Lagrangian. As it stands, the action in (|3.1q ) is given for an arbitrary choice of L and 
R consistent with the constraints. There is, of course, an enormous amount of redundant 
information contained in this description, since many L's and i?'s are equivalent to each 
other through a change of coordinates. To obtain an action which only depends on the 
truly physical variables p, f we make a specific choice for L and R, ie. choose a gauge. In 
so doing, we must respect the condition 



R'(f + e) - R'(r - e) 



1 

R 



p2 _j_ m 2£2 



which constrains the form of R' arbitrarily near the shell. Suppose we choose R for all 
r > f ; then R'(r — e) is fixed by the constraint, but we can still choose R for r < f — e, 
in other words, away from the shell. We will let R' K denote the value of R! close to the 
shell but far enough away such that R is still freely specifiable. We employ the analogous 
definition for R! > , except in this case we are free to choose R'y — R'(f + e). 

In terms of this notation the time derivative of S is 

. dS 



dt 



rRL 



-RRlog 



(R^/L) 2 - 1 + 2M/R - JiR^/L) 2 - 1 + 2M+/R 



R'{f - e)/L - J(R'(f - e)/L) 2 - 1 + 2M/R 



R'JL - J {R'JL) 2 - 1 + 2M/R 



+ / dr[ir R R + ir L L]+ dr[ir R R + n L L] - M+. (3.19) 

J r min J r+e 

At this point we will, for simplicity, specialize to a massless particle (m = 0) and define 
7] = ± = sgn(p). Then the constraints (|3.12|) read 



R'(f - e) = R'(f + e) + 



rjp 
R 



\R'(f -e)/L) 2 -1 + 2M/R 
These relations can be inserted into (|3.19|) to yield 

L = 'rRL 



(R'(r + e)/L) 2 -l + 2M + /R+-^. 

LR 



(3.20) 



(R^/L) 2 - 1 + 2M/R - ^(R^/L) 2 - 1 + 2M+/R 



-r/RR log 



R'JL - riJiR^/L) 2 - 1 + 2M+/R 



+ 



R'JL - 77 y (R^/L) 2 - 1 + 2M/R 
dr [ttrR + n L L] + / dr[n R R + n L L] - M 4 



(3.21) 



9 



Now we can use the freedom to choose a gauge to make ( p. 21 ) appear as simple as possible. 
It is clearly advantageous to choose L and R to be time independent, so tcrR + n L L = 0. 
Also, having R' = L simplifies the expressions further. Finally, it is crucial that the 
metric be free of coordinate singularities. A gauge which conveniently accommodates 
these features is 

L = 1 ; R = r 

The Schwarzschild geometry in this gauge is reviewed in the Appendix. It is considered 
in more depth in ||. 

The L = 1, R = r gauge reduces the Lagrangian to 



•[V2Mf 



2M + f] — t]rr log 



Mr 



•r - r]V2M 

where M + is now found from the constraints ( |3.20|) to be related to p by 

M.-M 



(3.22) 



P 



r] - ^2M + /r 

The canonical momentum conjugate to f obtained from |3.22| is 



(3.23) 



dL r— 

p c = — = V2Mf 
dr 



2M + f — rjf log 



, (3.24) 



(3.25) 



which identifies M + as the Hamiltonian. We should point out that M + is the Hamiltonian 
only for a restricted set of gauges. If we look back at Q3.21 ) we see that the terms hrR+ixlL 
will in general contribute to the Hamiltonian. 



in terms of which we write the action in canonical form as 



S= dt\pj-M+ 



4. Quantization 

In this section we discuss the quantization of the effective action (|3.25|). First, it is 



convenient to rewrite the action in a form which explicitly separates out the contribution 
from the particle. We write 

M + = M-p t 
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so 

S= ldt\p c r + p t ) (4.1) 



and the same substitution is understood to be made in (|3.24|) . We have omitted a term, 
J dtM, which simply contributes an overall constant to our formulas. In order to place 
our results in perspective, it is useful to step back and consider the analogous expressions 
in flat space. Our results are an extension of 



p = ±\fPt 2 — ra 2 (4.2) 

S = Jdt[pr+p t }. (4.3) 

Indeed, the G — > limit of ( |3.24| ), ( p.25| ) yields precisely these expressions (with 
m — 0). To quantize, one is tempted to insert the substitutions p — > —i-§^, Pt —> — 2J7 
into ( |4.2j ), so as to satisfy the canonical commutation relations. This results in a rather 
unwieldy, nonlocal differential equation. In this trivial case we know, of course, that the 
correct description of the particle is obtained by demanding locality and squaring both 
sides of ( [4.2| ) before substituting p and p t . So for this example it is straightforward to 



move from the point particle description to the field theory description, i.e. the Klein- 
Gordon equation. Now, returning to ( |3.24j ) we are again met with the question of how 
to implement the substitutions p — > —id. In this case the difficulty is more severe; we no 
longer have locality as a guiding criterion instructing us how to manipulate ( |3.24| ) before 
turning the p's into differential operators. This is because we expect the effective action 



(4.1) to be nonlocal on physical grounds, as it was obtained by including the gravitational 



field of the shell. 

There is, however, a class of solutions to the field equations for which this ambiguity is 
irrelevant to leading order, and which is sufficient to determine the late-time radiation from 
a black hole. These are the short-wavelength solutions, which are accurately described by 
the geometrical optics, or WKB, approximation. Writing these solutions as 

(j)(t,r) = e iS{t ' r) , 

the condition determining the validity of the WKB approximation is that 

|asi » \o 2 s\ 1 / 2 , id^i 1 / 3 ... 

and that the geometry is slowly varying compared to S. In this regime, derivatives acting 
on <p{t, r) simply bring down powers of dS, so we can make the replacements 

dS OS 

Pc "Tj ; Pt "777" 

or ot 
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and obtain a Hamilton- Jacobi equation for S. Furthermore, it is well known that the 
solution of the Hamilton- Jacobi equation is just the classical action. So, if r(t) is a 
solution of the equations of motion found by extremizing (|4.1|) , then 



where 



S(t,r(t)) = S(0,r(0)) + / dt p c (f(t))f(t) + p t 

Jo L 



dS 

Pc(0,r) = (0,f) . 



(4.4) 



(4.5) 



Since the Lagrangian in (|4.1|) has no explicit time dependence, the Hamiltonian p t is 
conserved. Using this fact, it is easy to verify that the trajectories, r(t), which extremize 
(|4.1| ) are simply the null geodesies of the metric 



ds 2 = -dt 2 + \dr + 



I2M, 



dt 



(4.6) 



From ( |A.4j ) the geodesies are: 



ingoing: t + f(t) + 2J2M + r(t) + 4M + log [Jf(t) + J2M. 



f (0) + 2 J2M+f(0) + 4M+ log [Jf(0) + J2M. 



outgoing: t - f(t) - 2y / 2M + f(t) - 4M+ log [y'f(t) - yj2M. 



= -f(0) - 2^/2M + f(0) - 4M+ log [yjf(0) - yj2M+]. (4.7) 
M + , in turn, is determined by the initial condition S(0,r) according to (|3.24 ) and ([4.5|): 



dS 



ingoing: —(0, f(0)) = J2Mr(0) - J2M+r(0) + r(0) log 
ar v v 

dS 



outgoing: — (0,f(0)) = J2Mr(0) - J2M + r(0) - f(0) log 



^/H0) + V2M 

^Jm-V2w; 



f(0) - v/2M 



(4.8) 



Finally, we can use this value of M + to determine p c (t): 



ingoing: p c (t) = yJ2Mr(t) - yj2M + f(t) + f(t) log 



outgoing: p c (t) = ^2Mr{t) - ^2M + r(t) - f(t) log 




(4.9) 
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These formulas are sufficient to compute S(t,r) given S(0,r). 

As will be discussed in the next section, the relevant solutions needed to describe the 
state of the field following black hole formation are those with the initial condition 



S(0,r) = kr k>0 (4.10) 

near the horizon. Here, k must be large 1/M) if the solution is to be accurately 
described by the WKB approximation. In fact, the relevant fc's needed to calculate the 
radiation from the hole at late times become arbitrarily large, due to the ever increasing 
redshift experienced by the emitted quanta as they escape to infinity. We also show in 
the next section that to compute the emission probability of a quantum of frequency u, 
we are required to find the solution for all times in the region between r = 2M and 
r = 2(M + uj). That said, we turn to the calculation of S(t,r) in this region, and with 
the initial condition ( |4.1U| ). The solutions are determined from fl4.4| ), fl4.7D -( fO| ). Because 
of the large redshift, we only need to keep those terms in these relations which become 
singular near the horizon. We then have for the outgoing solutions: 



S{t,r) = fcr(0) 



f(0) 



dr f log 



Vr- y/2~M+ 



(M+ - M)t 



t - 4M+ log [y/r - \J2M+ 
k = — f (0) log 



4M+ log [Jr(0) - J2M- 



f(0) - 



yr(0) - V2M 



(4.11) 
(4.12) 
(4.13) 



To complete the calculation, we need to invert (|4.12 ) and ( [4.1 3| ) to find M + and f(0) 
in terms of t and r, and then insert these expressions into ( |4.11| ). One finds that to next 
to leading order, 



2M+ = V2M + (Vr - V2M) 



yr(0) = V2M + {Vr- V2M) 



where 



M' M • v / 2M( v / r - V2M) 



I _)_ ( e k/2M' _ l^ e -t/4M" 

e (k/2M'-t/4M') 
I _|_ ^ e k/2M' _ lj e -t/4M' 

(k/2M-t/4M) 



(4.14) 



(4.15) 



1 _|_ e (k/2M-t/4M) ■ 

Plugging these relations into (|4.11|) and keeping only those terms which contribute to the 
late-time radiation, one finds after some tedious algebra, 



S(t,r) = -(2M 2 -r 2 /2)log 



1 + e 



(k/2M'-t/AM') 



(4.16) 
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5. Results 



We will now discuss the application of these results to the problem of black hole 
radiance. We begin by recalling some general features of the quantization of a scalar field 
in the presence of a black hole ||. The quantization proceeds by expanding the field 
operator in a complete set of solutions to the wave equation, 



4>{t,r) 



dk 



a k fk(t,r) + alf£(t,r) 



(5.1) 



There are, however, two inequivalent sets of modes which need to be considered: those 
which are natural from the standpoint of an observer making measurements far from the 
black hole, and those which are natural from the standpoint of an observer freely falling 
through the horizon subsequent to the collapse of the infalling matter. The appropriate 
modes for the observer at infinity are those which are positive frequency with respect to 
the Killing time, t. Writing these modes as 



u k (r)e 



(p(t, r) reads 



cj)(t,r) = J dk 
These modes are singular at the horizon 

du k 



a k u k (r)e-^ kt + 4?4(r)e ia;fei 



(5.2) 



dr 



oo as r 



2M. 



Symptomatic of this is that the freely falling observer would measure an infinite energy- 
momentum density in the corresponding vacuum state, 



(O t |TjO t ) -> oo 



as 



2M 



where d k \0 t ) = 0. However, we do not expect this to be the state resulting from collapse, 
since the freely falling observer is not expected to encounter any pathologies in crossing 
the horizon, where the local geometry is entirely nonsingular for a large black hole. To 
describe the state resulting from collapse, it is more appropriate to use modes which 
extend smoothly through the horizon, and which are positive frequency with respect to 
the freely falling observer. Denoting a complete set of such modes by v k (t,r), we write 



<f){t,r) dk \b k v k (t,r) + i>lv* k (t,r) 



(5.3) 



Then, the state determined by 



b k \0 v ) = 
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results in a non-singular energy-momentum density at the horizon, and so is a viable 
candidate. The operators d k and b k are related by the Bogoliubov coefficients, 



dk' 



otkk'bk' + Pkk'bk 



(5.4) 



where 



1 



OCkk' 



0, 



kk> 



2nu k (r) 
1 

2nu k (r) 



oo 

-oo 
oo 

-oo 



dt e iuJkt v k/ (t,r) 
dte luJkt v* k ,(t,r) 



The average number of particles in the mode Uk(r)e tuJkt is 



N k = (0 v \ala k \0 v ) = dk' \p, 



>kk' 



(5.5) 



(5.6) 



If we treat the black hole as radiating for an infinite amount of time, then Nk will be 
infinite. To obtain the rate of emission we can place the hole in a large box and use the 
density of states duj/2^ for outgoing particles. Further, if \ot kk > / j3 kk i\ is independent of k', 
as will be shown to be the case, we can use the completeness relation 



dk' 



\otkk'\ 2 - \P, 



kk' 



(5.7) 



to obtain for the flux of outgoing particles with frequencies between u k and u k + du k 

du k 1 



2vr \a kk > / j3 kk >\ 2 - 1' 



(5. 



This gives the flux of outgoing particles near the horizon. As the particles travel outwards, 
some fraction, 1 — T(u), of them will be reflected back into the hole by the spacetime 
curvature. Thus, for the flux seen at infinity we write 



Foo{uj k ) 



duoi 



r(a> fe ) 



27T \a kk i/(3 kk i 



(5.9) 



Next, we consider the issue of determining the modes v k (t,r). As stated above, we 
require these modes to be nonsingular at the horizon. Since the metric near the horizon is 
a smooth function of t and r, a set of such modes can be defined by taking their behaviour 
on a constant time surface, say t — 0, to be 



v k (0,r 



« e lkr as r 



2M. 



This is, of course, the initial condition given in fl4.10| ). Now, the integrals in ( |5.5|) deter- 
mining the Bogoliubov coefficients depend on the values of v k (t,r) at constant r. Since 
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v k is evaluated in the WKB approximation, the highest accuracy will be obtained when r 
is as close to the horizon as possible, since that is where v^s wavelength is short. On the 
other hand, in calculating the emission of a particle of energy u k , we cannot take r to be 
less than 2(M + u k ), since the solution Uk(r)e~ VJkt cannot be extended past that point. 
Therefore, we calculate the integrals with r = 2(M + u k ). 

The results of the previous section give us an explicit expression for v k . From ( 4.16Q , 

V k (t, 2(M + U k )) = e iS <tX M+U *» = e i(4Afa, fc+2w »)log[l+e(*/^'-*/^')] ( g 1Q) 

where M' is 

e (k/2M-t/4M) e (k/2M-t/4M) 



M' = M + V2M(J2(M + cu k ) - « M + cu k 



I _|_ e (k/2M-t/AM) ' K I _|_ e (k/2M-t/4M) ' 

(5.11) 

Then, the integrals are, 

rfte io; fc t e ±i(4A^ fc +2^)log[l+e( fc / 2M, -'/ 4M ')] 5 /g -^x 

the upper sign corresponding to a^', and the lower to (3 kk r. We can compute the integrals 
using the saddle point approximation. It is readily seen that for the upper sign, the saddle 
point is reached when 

e (k/2M'-t/AM') _^ ^ 

so t is on the real axis. For the lower sign, the saddle point is 

e (k/2M'-t/4M') _ -1/2, 

which, to zeroth order in u k , gives 

t = 4inM + real 

and to first order in uj k , gives 

t = 4in(M - u k ) + real. 

Inserting these values of the saddle point into the integrands gives for the Bogoliubov 
coefficients, 

Otkk' 



Pkk 

The flux of radiation from the black hole is given by 



e 4 " (M -" fe K (5.13) 
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There is an alternative way of viewing the saddle point calculation, which provides 
additional insight into the physical origin of the radiation. Let us rewrite the integral 
(151%) as 

/• fYl 

dte iu, k t±iS(t,2(M+u k )) _ ( 515 ) 



The saddle point is given by that value of t for which the derivative of the expression in 
the exponent vanishes: 

dS 

u k ±—{t,2{M + uj k )) = 0. 
But dS/dt is just the negative of the Hamiltonian, 

Pt = M-M+ 

so the saddle point equation becomes 

M + = M±u k . 

To find the corresponding values of t, we insert this relation into ( 4.12 ) and ( f4.13|) : 



dS_ 
~dt 



t = A{M± u k ) log 



k = — f (0) log 



2(M + u k ) 



2(M + oo k ) 



y/f{G)-y/2(M±U k ) 



f(0)- J2(M±u k ) 



f(o) - \FM 



(5.16) 



(5.17) 



where we have written f = 2{M + uj k ) + e to make explicit that f must lie outside the 
point where the solutions u k (r) break down. We desire to solve for t as k — > oo. For the 



upper choice of sign, we find from ( |5.17|) that 

V^(0) = y/2(M + u k ) + 0(e" fc / 2M ) , 
which, from ( |5.16| ), then shows that the corresponding value of t is purely real. 



For the lower choice of sign we have, 



y/m = y/2(M + u k ) - 0(e- fc / 2M ) 



Continuing t into the upper half plane, we find from (|5.16|) that 

t = Ain(M - uo k ) + real. 
These results of course agree with our previous findings. 
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The preceding derivation invites us to interpret the radiation as being due to negative 
energy particles propagating in imaginary time. The particles originate from just inside 
the horizon, and cross to the outside in an imaginary time interval 47r(M — u^). This, 
perhaps, helps clarify the analogy between black hole radiance and pair production in 
an electric field, which, in an instanton approach flOU , is also calculated by considering 
particle trajectories in imaginary time. 



6. Comments 



1. Let us return to the question of thermality. One might have guessed that the 
correct exponential suppression factor could be the Boltzmann factor for nominal tem- 
perature corresponding to the mass of the hole before the radiation, after the radiation, 
or somewhere in between. Thus one might have guessed that the exponential suppression 
of the radiance could take the form e _w / Tbcforc , e~ u ' Ta&er , or something in between. Our 
result, to lowest order, corresponds to the nominal temperature for emission being equal 
to T a f teT . 

2. Our methods clearly generalize in a completely straightforward manner to other 
forms of black holes (e. g. Reissner- Nordstrom or dilaton holes) and to charged or massive 
scalar fields. We have undertaken extensive calculations along these lines, whose results 
will be reported elsewhere. It should also be possible to consider emission into higher 
partial waves, though this involves some new issues in treating the collective modes for 
rotations of a non-symmetric hole, that we have not yet investigated seriously. Similarly, 
it would be interesting to consider emission with recoil. 

3. One can also consider the case of two mutually gravitating particles. This involves 
geometries with two shell discontinuities. By analyzing this problem, we expect to be 
able to address the question whether there are non-trivial correlations in the radiation, or 



between incoming particles and subsequent radiation. There are however claims [11], [12 
that semiclassical methods become internally inconsistent, or at least highly singular, in 
the latter case. 
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Appendix: Schwarzschild Geometry with L = 1,R = r. 



In this gauge the line element for the Schwarzschild geometry reads 

ds 2 = -dt 2 + (dr + \—dt) 2 + r 2 {dB 2 + sinfl 2 ^ 2 ). 



(A.l) 



These coordinates are related to Schwarzschild coordinates, 

ds 2 = -(1 - — )dt 2 s + dr 2M + r 2 (d6 2 + sm6 2 d(f) 2 ) 



1 - 



(A.2) 



by a change of time slicing, 

t s = t- 2v / 2Mr - 2M log 



(A.3) 



In contrast to the surfaces of constant t s , the constant t surfaces pass smoothly through 
the horizon and extend to the future singularity free of coordinate singularities. 

In terms of r and t, the radially ingoing and outgoing null geodesies are given by 
ingoing: t + r — 2V2 Mr + 4M log[Vr + V2M] = v = constant 
outgoing: t-r- 2v / 2Mr - AM log [y/r - V2M] = u= constant . (A.4) 
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